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ABSTRAGT 

A tidal radius is a distance from a satellite orbiting in a host potential beyond which its material is 
stripped by the tidal force. We derive a revised expression for the tidal radius of a rotating satellite 
which properly takes into account the possibility of prograde and retrograde orbits of stars. Besides 
the eccentricity of the satellite orbit, the tidal radius depends also on the ratio of the satellite internal 
angular velocity to the orbital angular velocity. We compare our formula to the results of two A^-body 
simulations of dwarf galaxies orbiting a Milky Way-like host on a prograde and retrograde orbit. The 
tidal radius for the retrograde case is larger than for the prograde. We introduce a kinematic radius 
separating stars still orbiting the dwarf galaxy from those already stripped and following the potential 
of the host galaxy. We find that the tidal radius matches very well the kinematic radius. Our results 
provide a connection between the formalism of the tidal radius derivation and the theory of resonant 
stripping. 

Subject headings: galaxies: dwarf — galaxies: interactions — galaxies: kinematics and dynamics — 
galaxies: structure 


1. INTRODUGTION 

The tidal radius rt is a theoretical boundary of a satel¬ 
lite object (e.g. a star cluster, a galaxy) beyond which 
its material is stripped by the tidal forces originating 
from the potential of a larger host object (e.g. another 
galaxy, a galaxy cluster). I t was proposed for the first 
time by Ivon Hoerneil (|I957f ) in the context of globular 
clusters orbiting around the Milky Way. A strict theoret¬ 
ical definition exists only for satellites on circular orbits, 
for which the tidal radius is identical to the position o f 
L1/L2 Lagr ange points fsee iBinnev fc Tremaind 1200811 . 
iKind (jl962t l argued that for an eccentric orbit of a satel¬ 
lite, a star located at the instantaneous tidal radius is 
stationary in a rotating reference frame and is attracted 
neither towards the satellite nor towards the host. He 
also argued that the satellites are truncated during peri- 
center passages to the size indicated by the pericentric 
tidal radius. 

However, _ alre ady iHeno 2 (liMill and 

iKeenan fc InnanenI (jl975D noticed that retrograde 
orbits in the restricted three-body problem are stable up 
to much l arger distances tha n prograde ones, e xceeding 
the iKinei (|1962ll tidal radius. iRead et all (1200611 derived 
an expression for rt that takes into account three basic 
types of orbits: prograde, radial and retrograde. Indeed, 
the tidal radius for retrograde orbits turned out to be 
larger than for radial ones, which in turn was larger 
than for prograde orbits. 

Interactions of galaxies with the potential of larger ob¬ 
jects have proved capable of transforming their morphol¬ 
ogy and kinematics on various scales: th ose of dwarf 
galax ie s orbiting the Milk y Way- like galaxy ()Maver et al.l 
120011: iKazantzidis et al.l 1201111 . g alaxies in groups 
(IVillalobos et al.ll2012ll a nd clusters (iMastropietro et ^ 
1^051 : iBialas et al.l 1201511 . One of the factors influenc¬ 
ing the final outcome of such a process is the initial 
incli nation of the galaxy disk with respect to its or¬ 
bit (jVillalobos et al.l [20121 : iLokas et al.l 120151 1. The re¬ 


sult of an encounter of equal-mass disk galaxies also de¬ 
pends on the direction o f their rotation ()HolmberBlll941l : 
iToomre fc Toom^ll972l l. 

The tidal r adius is irnportan t for studying glob- 
ular clusters (iWebb et all 1201311 and dwarf galaxies 
((Lokas et al.ll2013ll as it can help disentangle their bound 
parts from the mate rial which has been alre ady lost 
and forms tidal tails (iKupper et al.l[2M^ 1201211 . Semi- 
analytical models of galaxy formation also rely on com¬ 
parisons of the tidal radius to t he size of the galaxy 
to describe timescales of stripping (iHenrioues fc ThomasI 
IMotiGhang et al.ll20H . 

In this paper we derive a new expression for the tidal 
radius, which properly takes into account both the eccen¬ 
tric orbit of the satellite and the orientation of its internal 
rotation. We apply this formula to iV-body simulations 
of a disky dwarf galaxy orbiting a Milky Way-like host 
and compare its predictions to features detected in mass 
and velocity distributions of the dwarf. 

2. DERIVATION 

Let a dwarf galaxy with a spherically symmetric mass 
distribution m{x) orbit around a host galaxy with a 
spherically symmetric mass distribution M{x). In a ref¬ 
erence frame centered on the dwarf the acceleration of a 
star belonging to the dwarf, located at Xs, is given by 

= -Gm(|a:s|)-r^ - GM{\xs - a;h|)-4^—^ 

|Xs| \Xq X]q| 

-GMda^hl)^, (1) 

\XYi\ 

where Xh is the position of the host and G is the grav¬ 
itational constant. The first term is the gravitational 
attraction by the dwarf galaxy, the second one is the at¬ 
traction by the host galaxy and the last one is the inertial 
force arising because the reference frame is not inertial, 
as the dwarf is falling onto the host. The second term 
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can be expanded in the series for |a;s|/|a;h| 1, yielding 

|^h| 

( 2 ) 

where p(|a;h|) = dlogM/dlogx|a;=l^j^l is the logarithmic 
derivative of M{x) calculated at x = |a;h| and the dot 
denotes the scalar product. 

The second and third terms together are commonly 
referred to as the tidal force. It squeezes the dwarf in 
the plane perpendicular to the direction to the host and 
elongates it in the direction towards and away from the 
host. In such an environment an initially disky dwarf 
is lik ely to form a tidally induced bar (e.g. iLokas et al.l 

[Ml . 

Let us change the reference frame to one rotating with 
a variable angular velocity O. New terms, the inertial 
forces, appear: the Euler force, the Coriolis force and 
the centrifugal force. The acceleration of the star now 
reads 

Xs = -Gm(|a;s|)-^ - GM(|a;h|)-r-^ 

|®s| l^h| 

+ GM(|xh|)[3-p(|a;h|)] ^h)^h _ ^ aSg 

- 20 X is - ri X (11 X cCs) + G [(|a;s|/|a;h|)^] , (3) 

where the cross denotes the vector product. 

In order to proceed we need to make an assumption 
about the velocity of the star. Let us assume the star 
follows a circular orbit around the dwarf with angular 
velocity Og, so in the non-rotating reference frame cen¬ 
tered on the dwarf its velocity is 

is ~ ^ ^s- (4) 

In the rotating frame we have to subtract the velocity of 
the frame itself, therefore 

is = (Us — H) X Xs- (5) 

This expres sion is substantiall y different from what was 
assumed bv IRead et al.l (|2006ll and we discuss this issue 
below. After substitution of the velocity we obtain 

is = -GmdaJsD-j-^ - GM(|a;h|)-j-^ 

-b GM(|xh|)[3-p(|a;h|)] ^h)^h _ ^ a^s 

- 211 X (Us X Xs) + n X {n X Xs) + O [(|a:s|/|a:h|)^] , 

( 6 ) 

which is the final formula for the acceleration of the star 
on a circular orbit around the dwarf galaxy, calculated in 
the frame rotating with a variable angular velocity. We 
note that up to this point vectors Us and U are arbitrary. 

We now derive the tidal radius for stars whose orbits lie 
in the same plane as the orbit of the dwarf. We choose U 
so that the unit vector ih = a;h/|a:h|, pointing from the 
dwarf towards the host, is constant in time. Such an U is 
equal to the instantaneous orbital angular velocity vector 


of the dwarf. In such a setup H || Us and U T Xs- Using 
the identity A x {B x C) = [A ■ C)B — (A • B)C and 
these properties we can rewrite double vector products 
obtaining 

is = -Gm(|xs|)^-GM(|.H|)^ 

-b GM(|a;h|)[3 -_p(|a:h|)]^^Y"^TT^ - n x Xs 

\Xh\ 

+ 2(11 • Us)a;s - \ fl\‘^Xs + O [(|a;s|/|£Ch|)^] • (7) 

Stars are mainly stripped when they are at the smallest 
or largest distance from the host and leave the vicinity 
of the dwarf through the LI and L2 Lagrange points. To 
find the component of the acceleration along the direc¬ 
tion towards the host assuming Xg || a^h we multiply the 
above equation by Xg = a;s/|a;s| to get 


-L X 

Xg- Xs = -Gm{xs)— + GM{xh)[2 -p(a;h)]-|- 

-b ^2^ - Xs -b G [(Xs/Xh)^] , (8) 

where we denoted magnitudes of each vector as |a| = a, 
except for Ug which obeys the relation U • Ug = UUs, 
i.e. its sign depends on its direction with respect to the 
direction of U. For stars on prograde orbits with respect 
to the orbital velocity of the dwarf Ug > 0, while for stars 
on retrograde orbits Ug < 0. 

We can parametrize the orbital angular velocity by the 
relation to the circular orbital angular velocity of the host 
Ucirc(a:h) Q-S 

q 2 _ ^circ(^h) _ GM(xh) , . 

A(xh) a:^A(xh) ' 


The parameter A(xh) depends on the orbit of the dwarf 
and details of the host potential. For a circular orbit A = 
1, for a point-mass host A(xh) = Xh/[a(l — e^)], where a is 
the major semi-axis and e is the eccentricity of the orbit. 
One may wonder if we can further assume something 
about Us, for example that it is equal to the circular 

1 /2 

velocity in the dwarf, i.e. Ug^circ = ± [Gm(xs)/Xs] 
However, such an assumption leads to incorrect results, 
as we demonstrate below. 

The tidal radius rt is commonly defined as the distance 
from the center of the dwarf where there is neither ac- 
celeration t owards the dwarf nor towards the host (e.g. 
[Kind Il962t [Read et all 1200611 . In our case this can be 
written as Xg ■ ais|a;g=rt = 0. After applying this relation, 
substituting for and truncating the series we obtain 
the equation for rt 


-Gm{rt)\ + GM(xh)[2 -p(xh)]^ + 


^ GM(xh) 
x((A(xh) 



which can be rewritten in a more convenient form 

_ f [m(rt)/M(xh)]A(xh) 

V2a/U-l + [2-p(xh)]A(xh)y' ’ 

2Us/U > —[2 — p(xh)]A(xh) + 1. 


( 10 ) 

( 11 ) 


( 12 ) 

(13) 
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The above inequality has to be fulfilled, otherwise the 
tidal radius does not exist. In such a situation the accel¬ 
eration points towards the dwarf. One can interpret this 
as the tidal radius being infinite. However, we assumed 
that Tt <C hence such an extrapolation is not valid. 

If we assume the simplest setup, namely a point-mass 
host (so M{x\i) = M and p = 0) and dwarf (m(rt) = m), 
and a circular orbit of the dwarf (A = 1), we obtain 

/ 1 m \ 

’'■"-(wnTTMj • 

Hs/H>-l/2. (15) 


Assuming further that Hg/H = 1 we recover the ex¬ 
pression for the L1/L2 Lagrange points, sometimes 
also referred to as the Jacobi radius rj (see e.g. 
IBinnev fc Tremainell200^ : 



This is straightforward to understand: a star located at 
the L1/L2 point is stationary in the rotating reference 
frame, hence it rotates around the dwarf with the same 
angular velocity as the dwarf revolves around the host. 

From the above it is also clear why we cannot as¬ 
sume the star to be moving with angular velocity Hg = 
[Gm{xs)/. In such case the angular velocity at the 
true distance of the L1/L2 points would be Hg = 
thus significantly overestimated. Consequently, the tidal 
radius would be underestimated. Therefore, as ilg one 
has to use the angular velocity already perturbed by the 
host. 

It may seem that for a dwarf on a retrograde orbit 
(Hg < 0), which spins sufficiently fast, the tidal radius 
is infinite. However, as we have already argued, in this 
case our approximation breaks down. Moreover, in real 
objects |Og| decreases with the distance form the center, 
therefore there will exist a radius where the necessary 
inequality will be fulfilled. 

If we assume in equation (na) that the star is at the 
(instantaneous) Lagrange point (Hg/H = 1) we can re¬ 
cover two known expressions for the tidal radius. For the 
dwarf on a circular orbit (A = 1) aro und an extended host 
we rec over the formula 8.108 from iBinnev fc Tremaiii^ 
(j2008li . namely 


n 


Xh 


i(h) 


_3-p(a:h) M{xi,)_ 


1/3 


(17) 


For an elliptical orbit (A = a:h/[a(l —e^)]) around a point- 
mass ho st (jp = 0) we recover the limiting radius formula 
(7) from iKin^ (I1962D : 


rt = a^h 


m{rt)/M 


- 11/3 


(a/xh)(l - e^) -1- 2 


(18) 


We emphasize that the exceptionally simple formula rt = 
aih.peri {to/[M( 3 -I- e)]}^^^ is valid only at the pericenter 
of the satellite’s orbit (i.e. for Xh = aJh.peri = a(l — e)). 

Let us now comment on the diffe rences betw e en ou r 
derivation and the one presented in iRead et al.l (I2006D . 
Our formulae for t he star velocity ffl and m should 
correspond to what IRead et al.l (120061) assumed in their 


equation (3). In the non-rotating reference frame cen¬ 
tered on the dwarf their expression is equivalent to 

Ag = (Jig -)- H) X ajg. (19) 

This should be d irectly compared to our formula O- 
IRead et ^ (|2006[) also assume that ilg = aHg^circ, where 
f^s,circ is the circular angular velocity of the dwarf and 
a = 1, 0 or —1 accounts for, respectively, prograde, radial 
and retrograde or bits of the stars. Hence, according to 
IRead et aH (|2006ll . when the star is located at the tidal 
radius, it has an effective angular velocity of aHg^circ + fi- 
It means that the stars on prograde and radial orbits 
are accelerated by a constant value, equal to the orbital 
angular velocity of the dwarf. The stars on retrograde 
orbits are first decelerated and then accelerated, but then 
they orbit in a different direction than the body of the 
d warf. The r a dial o rbits, initially having Hg = 0, are 
in IRead et al.l (|2006D framework accelerated to Hg = H. 
Thus it is not surprising that the formula of IR ead et_aL 
(2006) reduced for the radial orbits to the one of iKing 
(I1962D . as equality of the angular velocities is precisely 
the assumption of the latter. 

From the tidal stirring simulations of dwarfs on pro¬ 
grade orbits it is known that the rotation of the dwarf 
dimin ishes, rather than in creases (see e.g. iLokas et al.l 
l2015|) . This explains why iLokas etld] ()2013D obtained 
good match of the tidal radius to A^-body simulations 
assuming radial orbits of the stars, whereas the dwarfs 
were in fact on mildly pr ograde orbit s . Simp ly, the stars 
did not rotate as fast as IRead et abl (|2006l) assumed in 
the prograde case. We discuss this issue further in the 
next sections. 

3. COMPARISON TO A^-BODY SIMULATIONS 

Any formula should be compared to a realistic situa¬ 
tion in order to check its validity and the assumptions 
made during the derivation. We shal l use two of t he N- 
body simulations studied in detail bv ILokas et al.l (120151) 
that followed the evolution of a dwarf galaxy around a 
Milky Way-like host. The live realiz ation of a Milky Way- 
like galaxy consisted of an NFW (iNavarro et al.l [19951) 
dark matter halo with a virial mass of 7.7 x 10^^ Mq and 
concentration c = 27 and an exponential disk of mass 
3.4 X 10^° Mq with a radial scale length 2.88 kpc and 
thickness 0.44 kpc. The dwarf galaxy had an NFW halo 
of 10® Mq and c = 20 and a disk of mass 2 x IO^Mq, 
the radial scale length 0.41 kpc and thickness 0.082 kpc. 
Each component of each galaxy contained 10® particles. 

The dwarf galaxy was initially placed at an apocenter 
of an eccentric orbit around the host galaxy with apoc¬ 
enter and pericenter distances of 120 and 25 kpc, respec¬ 
tively. The plane of the orbit was the same as the plane 
of the host’s disk. The only parameter which was varied 
between the simulations was the initial inclination of the 
disk of the dwarf with respect to the orbit. Here we will 
use only two configurations: t he exactly prograd e and 
the exactly retrograde, which in ILokas et al.l (I2015D were 
referred to as 10 and 1180, respectively. 

In principle, the tidal radius formula m applies to in¬ 
dividual particles in the dwarf, as each star has a different 
angular velocity Hg. In order to apply it to the dwarf as a 
whole, we measured the Hg(a:g) profile, i.e. the profile of 
the mean angular velocity component parallel to the or¬ 
bital angular velocity of the dwarf. Obviously, particles 
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Fig. 1.— Density profiles of the stellar component of the dwarf 
galaxy for the prograde model (red) and the retrograde model 
(blue), as well as dark matter density profile (black) at consecutive 
apocenters. 


in the dwarf are not on exactly circular orbits and have 
non-negligible radial velocities. As a result, an additional 
component of the Coriolis force appears. However, it is 
orthogonal to Xs and, consequently, does not have any 
impact on the force component along Xg. 

The right-hand side of equation (fT^ depends on rt, 
thus it is an implicit formula for the tidal radius, which 
can be solved numerically. In order to obey the assump¬ 
tions, from now on we will treat the disk of the host as a 
point mass. Furthermore, we calculated m{rt) assuming 
spherical symmetry of the dwarf galaxy mass distribu¬ 
tion. While these approximations introduce some error, 
we consider it negligible, as it only underestimates the 
gravitational force of the disks of both galaxies, which 
are at least an order of magnitude less massive than their 
dark matter haloes. The calculation of A(a:h) was carried 
out assuming that the dwarf is a point mass object or¬ 
biting a spherical host. 

In Figure [T] we show density profiles of the stellar com¬ 
ponent during consecutive apocenters for both models. 
In addition, we plot the profile of the dark matter com¬ 
ponent, which was very similar in both cases. In the case 
of the prograde model there is a clear break in the den¬ 
sity profiles where the main body of the dwarf ends and 
the tidal tails begin. On the other hand, the profiles of 
the retrograde model are rather featureless and close to 
a power law. 

For their mildly prograde models, iLokas et~aI1 (I2013f) 
compared the break radius rg to the theoretical tidal 
radius. Here, we were only able to locate the break in the 
pr ograde model a n d to measure it we used the method 
of I Johnston et al.l (|2002D . At each radius we calculated 
slopes of the density profile interior and exterior to that 
radius. Comparing the two slopes we searched for the 
innermost difference larger than the assumed threshold. 

Stripping of the stars should also influence their kine¬ 
matics. This is because in the first approximation the 
stars bound to the dwarf are mostly influenced by the 


Prograde model 


Retrograde model 
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Fig. 2.— Left column: The rotational velocity (red) and the 
radial velocity (blue) of the dwarf on the prograde orbit. Right 
column: The rotational velocity of the dwarf on the retrograde or¬ 
bit. Negative values indicate that the dwarf counter rotates with 
respect to its orbit. Rows of panels correspond to consecutive apoc¬ 
enters. 


dwarf potential, whereas the motion of the stripped ones 
is governed by the potential of the host. Therefore, we 
introduce the kinematic radius beyond which the kine¬ 
matics of the particles is no longer dominated by the 
potential of the dwarf. 

Particular realizations of this definition are different 
in both models. To illustrate this, in Figure [2] we plot 
profiles of the rotational velocity v^p = rHg of the dwarfs 
at the consecutive apocenters. The inner parts of the 
retrograde model counter-rotate with respect to the or¬ 
bital angular velocity, as indicated by the negative values 
of Vp. The maximum of the velocity curve slowly de¬ 
creases at subsequent apocenters due to mass loss. Par¬ 
ticles which were stripped and are already far away (more 
than 10 kpc) have positive angular velocity with respect 
to the center of the dwarf, as they orbit the host. Taking 
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Prograde model Retrograde model 



Fig. 3.— The angular velocity profile of stars (red), the circu¬ 
lar angular velocity cu rve (blue) and the effective angular velocity 
(green), as assumed bv IRead et al.l 1120091 . The left panel is for the 
prograde and the right one for the retrograde model. Negative val¬ 
ues in the retrograde model indicate counter-rotation with respect 
to the orbital angular velocity. Both panels show measurements 
from the simulations at the first pericenter passage, after 1.15 Gyr 
of evolution. 


this into account, we define the kinematic radius as the 
one where fig = 0. 

In the case of the prograde model, the angular veloc¬ 
ity in the inner parts decreases significantly as the dwarf 
is transformed and the random motions start to domi¬ 
nate (see e.g. Fig. 2 in lLokas et al.l[2M5l) . In the outer 
parts, v^p grows similarly as in the retrograde model. Un¬ 
fortunately, there is no such a striking feature as in the 
retrograde case. One can see a minimum in the outer 
parts, but it is very wide and shallow, if present at all, 
therefore it is not a suitable proxy for the kinematic ra¬ 
dius. Instead, we define rk as a radius where the radial 
velocity is equal to the rotational velocity, i.e. Vr = v^,. 
We note that there is a region in the profiles of Vr where 
Vr < 0, indicating that some particles are reaccreted onto 
the dwarf. 

Let us now assess the assumption of iRead et al.l (l2006f ) 
concerning the effective angular velocity of the stars in 
the dwarf. In Figure [3] we compare the profiles of the 
measured angular velocity of the dwarf fig, the circular 
angular velocity ns,circ and the effective a ngular velocity 
afls,circ+f^ assumed bv iRead et al.l (|2n06[) . The measure¬ 
ments were done at the first pericenter passage. Near the 
centers, the fig profile is shallower than ns,circ partially 
due to the former being measured in spherical b i ns. It 
is well visible that the assumption of iRead et al.l (l2006f l 
is not valid, as it significantly overestimates the angu¬ 
lar velocity, which increases the Coriolis force term and 
leads to a large underestimation of the tidal radius. The 
magnitude of fig is smaller than flg^circ even at the first 
pericenter, and later on the discrepancy is even larger. 
Thus, as we discussed previously, ns,circ cannot be used 
as a proxy for fig. 

Figure m illustrates the differences between the distri¬ 
butions of the stellar component in both models and their 
relation to the tidal radius calculated from our equation 
m- During the third pericenter passage, the dwarf on 
the prograde orbit is already significantly smaller than 
the one on the retrograde orbit. In the prograde orbit 
case, the material is also more effectively removed from 
the vicinity of the dwarf. The calculated tidal radii cor¬ 
respond to the actual sizes of the dwarfs. 

In Figure [5] we show tidal radii calculated from for¬ 
mula (HI and compare them to tidal radii calculated 
with different assumptions. Let ns recall that the Jacobi 


Prograde model 



Fig. 4.— Surface density maps of the stellar component of dwarf 
galaxies just after the third pericenter passage (5.5 Gyr). The green 
circles indicate the tidal radii calculated from formula (|12|l . The 
top panel is for the prograde and the bottom one for the retrograde 
model. 


formula (HU) is the most widely used expression for the 
tidal radius and corresponds to the distance of L1/L2 La¬ 
grange points, as if the dwarf was on a circular orbit of 
radius equal to its current distance from the host. In the 
instantaneous Lagrange point approximation we assume 
Dg/D = I and use our formula (fT31) . This is equivalent 
to the iKind (119621 ) approximation, only without the as¬ 
sumption of a point-mass host. 

The tidal radius for the retrograde model is always 
larger than for the prograde model. However, the dif¬ 
ference is smaller later on because the rotation of the 
prograde dwarf slows down with time. Small peaks in 
ct measurements around the pericenters are caused by 
the Coriolis force originating from the movement of the 
whole dwarf on such an elongated orbit that locally ex¬ 
ceeds the tidal force. This behavior might be different if 
we included higher order terms in the tidal force. 

In the case of the prograde model, the simpler approx- 
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Tidal radii 



t (Gyr) 

Prograde model 



t (Gyr) 

Retrograde model 



Fig. 5.— Top: Tidal radii calculated for the prograde (red) and 
the retrograde (blue) model. Middle: Tidal radius (red), Jacobi ra¬ 
dius (green) and instantaneous Lagrange point (orange) for the pro¬ 
grade model. Bottom: Tidal radius (blue), Jacobi radius (green) 
and instantaneous Lagrange point (orange) for the prograde model. 


imations give similar results. This is due to the fact that 
the ratio rig/n varies in a range from about 0.3 to 2, so 
the approximation of the instantaneous Lagrange point 
is not that far from reality. The Jacobi radius performs 
worse, as it assumes an instantaneous circular orbit of the 
satellite. However, in the case of the retrograde model 
the difference is huge and results from taking into ac¬ 
count the orientation of the dwarf’s rotation. 

In Figure |6] we compare our tidal radius with the kine¬ 
matic and break radii obtained from iV-body simulations 
for the prograde and retrograde model. The kinematic 
radius agrees remarkably well with the tidal radius. How¬ 
ever, the agreement is worse between the fourth and fifth 
pericenter. In the prograde case it is because the stellar 
component has already significantly evolved from the ini¬ 
tial disky configuration. In the retrograde case, the pro¬ 
file of the angular velocity in the outer parts is very flat 
and it does not cross zero up to a large radius, possibly 
due to numerical errors. On the other hand, rg is only 
a rough est i mate of the tidal radius. As discussed by 
iLokas et~an (|20I3D , it is because the density distribution 
needs time to respond to the external forces. 


Prograde model 



Retrograde model 



Fig. 6.— Top: Tidal radius (red), kinematic radius (green) and 
break radius (blue) for the prograde model. Bottom: Tidal radius 
(red) and kinematic radius (green) for the retrograde model. 


4. CONCLUSIONS 

In this work we studied the tidal radius of satellites 
orbiting in a host potential. We derived an improved 
formula for rt for stars revolving around the satellite in 
the same plane as the satellite orbits the host. Our for¬ 
mula (fT^ properly takes into account that the star might 
be on a prograde or retrograde orbit. We verified our for¬ 
mula against A^-body simulations of a disky dwarf galaxy 
on a co-planar orbit around a Milky Way-like host. The 
tidal radius corresponds to the kinematic radius, i.e. the 
radius beyond which kinematics of the stellar particles 
is dominated by the host potential and the particles no 
longer orbit the dwarf galaxy. 

We showed that the assumption made bv IRead et al.1 
ll200^ concerning the effective angular velocity of the 
stars in the dwarf is not fulfilled in the case of an elon¬ 
gated orbit of the dwarf. Unfortunately, no justification 
of th e used assumption is presented in the discussed pa¬ 
per. IRead et al.l (|2006[) claim that for short timescales 
they obtained excellent agreement between their formula 
and simulations. However, at their tidal radii one can 
notice only some depletion of stars, but certainly not 
a boundary of the system. Thus, in our opinion, their 
tidal radius is significantly underestimated. Moreover, 
they take into account particles on polar orbits, which 
are not well described by either their or our formalism. 

As already noticed by iLokas et~aI1 (I20I3D , the peri¬ 
centri c tidal radius solely, as calculated from e.g. iKind 
(|I962D simple formula, is not enough to describe the ex¬ 
tent of the satellite as it orbits around the host object. As 
both measured quantities, the break and the kinematic 
radius, show, the size of the satellite changes along the 
orbit. After the pericenter passage, t he dwarf ex pands 
again (see Figure [Bl). Therefore, the iKin i (fl^ con¬ 
jecture that the satellites are trimmed at the pericenter 
and then remain unchanged, is not valid. Hence, to de- 
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scribe properly the size of a satellite one should use more 
advanced prescription. 

Applying our formula (1121) to the observations may be 
difficult, as one needs to know the orbit of the sat ellite. 
Howev er, this is also needed if one wants to use iKind 
(jl962li expression. As an improvement with respect to 
the classical formula for the Jacobi radius (HSl) one can 
use formula (1141) , which holds under similar assumptions 
as the former (i.e. a circular orbit and a point-mass host), 
but takes into account the rotation of the satellite, which 
can be readily measured for many objects. In this simple 
approximation, as the orbital angular velocity one may 
assume the circular orbital angular velocity. 

Our derivation depends on two basic assumptions: the 
stars orbit the dwarf on circular orbits and they are 
stripped in the vicinity of the line connecting the dwarf 
and the host (i.e. near L1/L2 Lagrange points). While 
the latter seems to be fulfilled based on the analysis o f 
tidal tails (|Klimentowski et ^120091: iLokas et al.lfeoiSl), 
the fo rmer is quite disputable. For example. iSzebeheld 
(|1967ll showed that in the restricted three-body prob¬ 
lem prograde orbits of the smallest body get complicated 
when particles approach the Lagrange points. On the 
other hand, retrograde orbits are fairly circular. Thus, 
we expect that our calculations are more robust in the 
retrograde case. 

Using an i mprov e d i mpulse approximation, 
iD’Onghia et al.l (j2009l 120101) proposed that strip¬ 
ping of the satellites is of resonant origin. Namely, 
for prograde encounters particles have largest velocity 
gains if their a^es = Hg/H ^ 1 during the pericenter 


passage. ILokas et'^D (|2015D calculated radii r^es at 
which Pres = 1 for the prograde model. In the vicinity 
of the pericenter, where most of the stripping occurs, 
Tj-es < i-O. resonant radius is smaller than the tidal 
radius. 

In view of the above, we may interpret the mechanism 
of stripping in the case of coplanar encounters as follows. 
In the prograde case the material in the outer parts of 
the disk is accelerated and pushed away by the resonant 
mechanism. When it travels beyond the tidal radius it 
becomes unbound from the dwarf and starts to follow the 
potential of the host. Some of this material is reaccreted 
by the dwarf, as its tidal radius increases when the dwarf 
recedes from the pericenter. In the retrograde case stellar 
pa rticles are accelerated much more weakly (see Fig. 4 
of ID’Onghia et al1l201(T[ ). so they are driven out of the 
dwarf, but they remain in its vicinity, as confirmed by 
their density profiles (Fig. [T]). 

Our analytic calculations are based on simplifying as¬ 
sumptions, which are not able to grasp other possible 
inclinations of the dwarf’s disk and the variety of par¬ 
ticle orbits. More investigation is needed to establish 
a proper connection between the stability of individual 
particle orbits and the general behavior of the density 
distribution and stripping. 


This work was supported by the Polish National Sci¬ 
ence Centre under grant 2013/10/A/ST9/00023. We ac¬ 
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